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We consider a model for superconductivity in a two-band superconductor, having an anisotropic 
electronic structure made of two partially overlapping bands with a first hole-like and a second 
electron-like fermi surface. In this pairing scenario, driven by the interplay between interband Vi,j 
and intraband Vj,i pairing terms, we have solved the two gap equations at the critical temperature 
T = T c and calculate T c and the chemical potential fi as a function of the number of carriers n for 
various values of pairing interactions, Vi,i, V?,2, and Vi,a. The results show the complexity of the 
physics of condensates with multiple order parameters with the chemical potential near band edges. 

I. INTRODUCTION 

In a metallic system where both a first band V\ and a second band V2, with no hybridization, cross the Fermi level 
the pairing is described by a matrix where the diagonal elements are the intraband Cooper pairing processes, namely, 
Vi,i and Vij,2 & n d the non-diagonal elements, namely, V\ 2 that describes the interband pairing process associated with 
the transfer of pairs between the bands [H3] . 

Today the open problem in condensed matter physics is the interband pairing where the standard BCS approxi- 
mation that considers the Fermi level far away form the energy of band edges or critical points breaks down. In fact 
where the chemical potential is tuned near an electronic topological transition the high temperature superconductivity 
has been associated with quantum interference effects, called shape resonances, in the interband pairing terms since 
1993 HE]. 

The shape resonance in the superconducting gaps is realized by material design of superconducting hetero-structures 
at atomic limit made of superlattices of superconducting layers with a relevant transversal hopping between layers 
[51 [7]. The material architecture and the superlattice misht strain [8 10 are key terms for fine tuning the system to 
quantum critical points [TT1 112| driven by electronic topological transitions. 

At the time of discovery of high Tc cuprates in 1986, these materials have been considered to be isotropic 3D 
perovskites but later it was recognized that all high temperature superconducting cuprates are made of a superlattice 
of atomic copper oxide layers intercalated by different spacers where the evolution of the complex electronic structure 
with chemical doping has remained mysterious for many years. The large majority of scientific community has 
interpreted the available unconventional experimental data in terms of an effective single band for high temperature 
superconducting models. Few authors have interpreted the data in terms of mainly two non-hybridized bands and 
two Fermi surfaces at the nodal and anti-nodal point of the band structure [5J [T3] ■ This key point for the physics 
of cuprates has been recently accepted by the large scientific community in 2010 driven by compelling experimental 
data [T3] . 

Later different practical realizations of superlattices of superconducting atomic layers intercalated by spacers tuned 
at a shape resonance [5J have been found in 2001 in diborides 15-30; and in 2008 in pnictides [3"TTl33| . 

In these hetero-structures many experiments have provided experimental tests for multigap superconductivity in the 
clean limit also in the presence of many lattice impurities introduced to change the chemical potential and structural 
phase separation from nano-scale to micron-scale [34H36] with the signatures of complexity of condensed matter near 
a quantum critical point. The key point for high temperature superconductivity is the presence of shape resonances in 
the interband pairing due to pair transfer between condensates with different topological symmetry due to quantum 
interference effects between pairing channel. This scenario has been discussed using simple electronic models [37J 13H] 

Superconducting pnictides [31] are made of stacks of iron fee atomic layers intercalated by spacers with a superlattice 
architecture similar to cuprates [7J that show multigap superconductivity [33j with multiple bands crossing the Fermi 
level [3J. 

While in diborides anomalous strong electron-phonon intraband interaction is quite relevant in two-band supercon- 
ductivity scenario, in pnictides the intraband coupling is so small that has been neglected [33]. In fact experimental 
findings show that the scenario where interband pairing is the dominant interaction, proposed by some of us in 
previous works 39 f [10], has now an experimental test. 
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Here we investigate the case where the relative strength between the intraband and interband coupling is changing 
that is relevant for material design of novel high Tc superconductors made of synthetic superlattices. In these systems 
the relevant region for high temperature superconductivity is where the chemical potential is tuned between the band 
edge of the second band and the electronic topological transition where the Fermi surface of the second band changes 
from a 3D topology to a 2D topology with a corrugated cylindrical shape [T5J HI] . Therefore we focus in a system 
with two similar overlapping bands such that the first hole-like (electron-like) Fermi surface coexists with a second 
electron-like (hole-like) Fermi surface. 

The proposed model has the characteristic feature that for the occupation number of 1 electron per spin shows a 
critical point where the hole Fermi surface is the same as the electron Fermi surface (like in pnictides [33 ) therefore 
tuning the chemical potential at this point there is a strong nesting between the two Fermi surfaces that determines a 
structural instability with a structural phase transition from tetragonal to orthorhombic phase |10j . A scale invariant 
complex structure could appear in the proximity of this phase transition in pnictides |10) as in cuprates |34j . 

We study a band structure that can be realized in superlattices of fee lattice like in iron based multilayer materials 
that have provided a particular realization of multiband superconductivity in the clean limit with a high critical 
temperature. 

In this paper we adopt the point of view that both the interband and intraband scattering terms are indeed very 
important to explain the response of the superconducting phase to the variation of the chemical potential. On the 
other side the case of pnictides is considered to be the case where intraband coupling can be neglected as we have 
proposed in the previous work before the discovery of pnictides [JD]- Now we focus our interest on the fact that 
it is important to consider different cases where both intra-band and interband pairing are relevant that will be 
relevant for novel synthetic superconducting superlattices. We consider the case where the interband pairing V\ \ is as 
important as the intraband contributions, namely, V 2<2l V\ %■ Here we report the study of the variations of the critical 
temperature and gap ratio as a function of the charge density crossing electronic phase transitions. 

This paper is organized as follows. In Section[ll]we present the four self-consistent equations: (1) two self-consistent 
equations (Eqs. ([6] - [7]) ) to be solved at T = T c , namely, the equation for the superconducting critical temperature 
T c coupled to the equation of the chemical potential /i at a given particle density n, and (2) three self-consistent 



equations for the gap parameters at T — and the total density (Eqs. (14 - 15 1). In Section III we present our 



numerical results. Finally we conclude in Section IV 



II. THE EFFECTIVE HAMILTONIAN AND THE BCS EQUATIONS 

The model Hamiltonian we use is the following 



k.k' 



at a} - a,-, .a ,-/ 



k.k' 



V 2 2 V bl 6 f - bp .b i (1) 



k,k' 

where 7 — ± =|, I and 

ex{k) = ei(fe)-/i (2) 

e 2 (fc) = e 2 (k)-n (3) 

are the three-dimensional tight-binding bands of bands 1 and 2. In previous papers [SHI 00], we have considered a 
Hamiltonian with V\ t % = V, and Vj.,1 = V 2 ,2 = 0, as the interband and intraband pairing potentials, respectively. 
Here, these conditions are relaxed. We consider the simple case of a cubic lattice 42 . The energy band dispersion is 
given by 

ei(fc) = ej — 2 1\ [cos(fc x ) + cos(fc y ) + S\ cos(fc 2 )] (4) 
e 2 {k) = e" - 2 £ 2 [cos(fcr) + cos(fcy) + s 2 cos(fc z )] . (5) 
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The Hamiltonian given in Eq. ([I]) is the same used in Ref. |39]s. 

Going through the same technical steps followed in Ref. [39] we obtain the following self-consistent equations 

1-n = F 1 {fi,T c ) 
1 = F 2 (^,T C , w c ) 



(0) 
(7) 



where 



F^,T C ) = Fi,i(/i,T c ) + F li2 (/J,T c ) 
F 2 (h,T c ,lo c ) = [v£ 2 -V 1A V 2 , 2 ] F 1 F 2 + V 1 , 1 F 1 + V 2 ,2F 2 



and 



FiA^T c ) 



dk 
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ex (A) 



" d 3 fc X2W 

e 2 (fc) 



2tt 3 



£i(fc) 
2/c s T c 

' £ 2 (fc) 
2k B T c 

tanh 
tanh 



gi(fc) 
2fc B T c 

g 2 (fc) 
2fc s T c 



(8) 
(9) 

(10) 

(11) 
(12) 

(13) 



where kg is the Boltzmann constant and Xi(fc) = 1 if | < w c and zero otherwise, i = 1,2 and w c is an energy 

cutoff (the Debye frequency) for the effective interaction. We solve the BCS equations at T = OK. These equations 
are: 



where 



G t = 

Ei{k) = 
A* = 
A? = 



Ai = AJGi 
A 2 = A^G 2 



77 £i(fc) £ 2 (fc) 

2tt 3 [^(fc) £ 2 (£) 



27T 3 
£i (fc)2 + (A:)2 

Vi,i Ai + Vl )2 A 2 
Vi, a Ax + V 2 . 2 A 2 



1,2 



1,2 



(14) 
(15) 
(16) 



(17) 

(18) 
(19) 



Results for the numerical solution of the self-consistent Eqs. (6|15 ) are presented in Section III 



III. SEARCHING FOR SUITABLE PAIRING PARAMETERS AND NUMERICAL RESULTS 

In Fig. fl] the two density of states Ni and N 2 , corresponding to the energy dispersion of Eqs. ([2]|3| are presented 
in panel (a) as a function of the electron number density per spin n for the band structure parameters chosen as 
ti = t 2 = 1 eV, a i = s 2 = 0.5, el — 0.0 eV and e 2 = 4.5 eV, i.e., two similar overlapping bands so that the first 
hole-like (electron-like) Fermi surface coexists with a second electron-like (hole-like) Fermi surface. 
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The region where the interband pairing is relevant is shown in panel (b) of Fig. [I] (for Vi,i = V 2i2 and Vi^/Vi.i = 
0.0,0.5,1.0,1.5 and 2.0) by an effective density tracer given by the ratio r(n) = V12 \/ Ni N 2 / A_ where A_ is the 
effective density of states given by the standard BCS theory as: 



1.13 lo c exp(A±) 



where 



A+ = 



2N X N 2 [Vi,iVi 



2,2 



JVi Via + N 2 V 22 ± N 



(20) 



AT = sj \Ni V M - N 2 V 2 . 2 f + 4 Aq N 2 



Finally, in panel (c) of Fig. [I] we present the self-consistent numerical solution of Eqs. (6|7) for chemical potential 
\x as a function of carrier density n which presents a quasi-linear dependence with n. 

Fig. [2] presents the self-consistent numerical solution of Eqs . (fl4 -16 1 for the energy gap ratio A2/A1 at T = OK 
in panel (a) and the self-consistent numerical solution of Eqs. (I6J|7| for T c x n at T = T c in panel (b). For these plots 
we have chosen Vi,i/t = V 2 . 2 jt = 0.50, w c /t = 1.00 and V h2 /t = 0.375, 0.750, 1.125, 1.500 and 1.875. 



IV. CONCLUSIONS 



We have provided a theoretical description for multiband superconductivity in novel synthetic superlattices of su- 
perconducting layers with cubic lattice structure that will be essential for material design of possible room temperature 
superconductors [6, 38!. We have investigated the two band superconductivity scenario for a two band superconduct- 
ing metal by an extended BCS mean field approach, considering the roles of the interband and intraband pairing. 
We have evaluated the critical temperature for different charge densities and chemical potentials, tuning the super- 
conducting system close to relevant energies for electronic topological transitions. The results indicate that both 
interband and intraband effective pairing are relevant in determining the properties of the superconducting phase and 
their sensibility to the variation of external parameters (i.e., density and chemical potential). Interestingly, for fixed 
intraband pairing, and several different interband pairing, the critical temperature as a function of the charge density 
displays characteristic resonances (and hence amplification of the critical temperature itself) when a critical point of 
the band dispersion is crossed by the chemical potential. Therefore, the systematic analysis of the critical temperature 
and the gaps can indicate a possible root to find new synthetic superconducting superlattices with high values of the 
transition temperature. Further work in progress investigating the isotope effect in this complex scenario. 
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iron pnictides and for quantum size effects in superlattices made of layers with cubic structure. The band structure for the 
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MgB2 hexagonal lattice (in 3-D) is given by 

1/2 



Ei,2 (*) = ±Vi ( A(fe) + B(k) 

+ 4 Va cos(fcito/2) cos(\/3fc y a/2) 

+ 2 V 3 cos(fc z c) (21) 

A(jfe) = [cos^a/v^) + 2 cos(fc !/ a/(2V3)) cos(fc^a/2)] 2 

B{k) + [2 sin(fcj,a/(2\/3)) cos(fea/2) - sin(fcs,a/V3)] 2 

where Vi.,2 is the hopping integral between atoms placed on the hexagonal lattice (second type of atoms placed inside the 
primitive cell). V3 is the effective hopping integral along the c-axis (or z— axis). The hexagonal lattice (Eq. ( 2 1 1 ) is left for 
a future calculation. 
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FIG. 1: Panel (a): the density of states Ni and N2 as a function of the carrier number n according to the chosen band 
parameters t\ = t% = 1 eV, si = S2 = 0.5, e? = 0.0 eV and e% = 4.5 eV; panel (b): the ratio r(n) = V12 V Ni N2 / A_ for 
= V2,2 and Vi^/Vi,! = 0.0,0.5, 1.0, 1.5 and 2.0 (see Eq. (201 for the definition of A_), useful as a tracer of the region where 
the interband pairing is relevant; panel (c): the self-consistent numerical solution of Eqs. ( 6jT I for chemical potential fi as a 
function of n at T — T c , which proved to be i ndependent of the choice of Vi,i, Vi,2 and u c and matches to the [s, x n 
relation at T = K obtained from Eqs. ( 14 - 16 1 (coloured online). 
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FIG. 2: For the given set of values Vi,i/t = V 2 , 2 /t = 0.50, uj c /t = 1.00 and Vi, 2 /t = 0.375,0.750,1.125,1.500 and 1.875 the 
panel (a) shows the self-consistent numerical solution of Eqs. ( |14| - |16[ l for the energy gap ratio A2/A1 at T = OK and panel 
(b) shows the self-consistent numerical solution of Eqs. (JsJ[tJ) for T c x n at T = T c . The solution of T c x n for Vip/t = 0.375 is 
not shown in panel (b) since the values for T c in this case are too close to zero (coloured online). 



